Revisiting Surface Diffusion in Random Deposition 
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An investigation of the effect of surface diffusion in random deposition model is made by analytical 
methods and reasoning. For any given site, the extent to which a particle can diffuse is decided by 
the morphology in the immediate neighbourhood of the site. An analytical expression is derived 
to calculate the probability of a particle at any chosen site to diffuse to a given length, from 
first principles. Using the method, the probabilities for different diffusion lengths are calculated 
and their dependence on system size and the number of deposited layers is studied. Numerical 
simulation of surface diffusion in random deposition model with varying extents of diffusion are 
performed and their results are interpreted in the light of the analytical calculations. Thus, a 
clearer understanding of the diffusion process and the effect of diffusion length on surface roughness 
is obtained. Systems with surface diffusion show nearly random deposition-like behaviour upto 
monolayer deposition. Their interface widths, in a logarithmic plot, are initially linear, as in random 
deposition. With increase in the number of layers, correlation effects between neighbouring columns 
become dominant. The interface deviates from its initial linear growth and eventually becomes 
saturated. An explanation for this behaviour is discussed and the point of departure from the linear 
form is estimated analytically. 

PACS numbers: 68.35.Fx, 68.35.Ct, 68.35.-p 
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I. INTRODUCTION 

Surface diffusion in random deposition model was in- 
troduced by Familjii to study the scaling properties of 
rough surfaces under the effect of surface diffusion in 
(1 + 1) dimensions. In the model, particles arc deposited 
from above onto a line of sites. The particle arriving at 
any randomly chosen site, diffuses to the nearest neigh- 
bouring site that is lower than the selected site. The 
resulting surface is smooth compared to the surface in 
random deposition. Diffusion to neighbouring sites in- 
troduces correlations among the neighbouring columns 
which eventually lead to saturation of the surface rough- 
ness. 

The roughness of the interface is quantitatively ex- 
pressed in terms of interface width W , which is a function 
of system size L and the growth time t. W{L,t) is de- 
fined by the root-mean-square (rms) fluctuation of the 
height of the interface. 
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Following Family- Vicsek^ scaling relation, W{L,t) scales 
initially as t^ for a fixed system size. It saturates to a 
time- independent value Wsat, beyond a mo del- dependent 
saturation time tsat- The saturated width scales with the 
system size L as Wsat ~ i"- 

It was observed by Family, through simulations in 
one dimension, that the scaling exponents obtained are 
quite different from those obtained for random deposi- 
tion (RD)^. An interesting point, reported by Family, 



was that the surface properties do not change with fur- 
ther increase in diffusion length. However, no explana- 
tion was put forward for this independence. 

Surface diffusion is found to be one of the most im- 
portant factors affecting surface morphology in several 
practical applications. For example, in molecular beam 
epitaxy (MBE), selective area metalorganic vapour phase 
epitaxy (SA-MOVPE) and chemical beam epitaxy of 
nanowires, the length of surface diffusion is significant in 
determining the growth rate and the mode of growth^. 
This diffusion length may be controlled by varying one 
or more experimental parameters such as the substrate 
temperature or introducing atomic or molecular hydro- 
gen under pressure^. An investigation is hence neces- 
sary to understand the role of surface diffusion in surface 
growth processes, and the nature of the dependence. 

In Family's model, the extent of diffusion is confined 
to nearest neighbours. However, at any stage of deposi- 
tion, there may be several sites with successively lower 
neighbours beyond the nearest. At such sites, a priori it 
seems, roughness may be further reduced if diffusion is 
not limited to nearest neighbour alone. 

In this present work, we study in an analytical manner 
and also through numerical simulation, the significance of 
diffusion length in determining interface roughness. Us- 
ing simple combinatorial argument, the probability with 
which a particle at any chosen site will diffuse over a cer- 
tain distance is derived, whence, for a given system size 
and number of deposited particles, the probability of dif- 
fusion over different distances (nearest neighbour, next 
nearest neighbour etc.) can be estimated. The diffusion 
length with the maximum likelihood will be the most 



dominant in determining the smoothness of the surface. 

To complement the analytical study, numerical sim- 
ulations of random deposition with surface diffusion are 
performed. A newly arriving particle is allowed to diffuse 
over a prescribed neighbourhood defined by a diffusion 
length Nd- Starting with N^ ~ 1, corresponding to Fam- 
ily's nearest neighbour diffusion, this diffusion length is 
varied over the values 2,3,4,6,.. .16 and its effect on the 
dynamics of growth is observed. We shall refer to this ex- 
tended Family's model as random deposition with surface 
diffusion (RDSD). The system sizes chosen, range from 
small sizes 24, 48, ■ • • 384, studied by Family, to larger 
sizes like 1024, 2048 and 3072. 

Our analytical study shows that the nearest neighbour 
diffusion probability is most dominant, thereby explain- 
ing Family's observation that diffusion to farther neigh- 
bours does not affect surface properties. The numerical 
simulations show a small change in surface properties as 
the diffusion length is increased from 1, 2 to 3. Beyond 
Nd = 3, the surface properties show no further change, 
in accordance with our analytical results. 

Wc address another related feature of RD and RDSD 
processes in the sub-monolayer region. For RD, the loga- 
rithmic plot of the interface width with time increases lin- 
early. With the introduction of surface diffusion, the plot 
is nearly collinear with that for RD in the sub-monolayer 
region. A deviation from this linear form is evident as 
a complete layer of deposition is approached, much be- 
fore the saturation of the interface width. This feature is 
discussed and a quantitative estimate of the time of this 
departure from the straight line is made. 



II. THE MODEL 

In Family's model of surface diffusion in random 
deposition^, a particle dropped in column i sticks to the 
top of the column i, i + 1 or i — 1 depending on which of 
the three columns has the smallest height. If the columns 
i, i + 1 and i — 1 have the same height, then the particle 
simply stays on top of column i and docs not diffuse. If 
the neighbouring columns i + 1 and i — 1 arc of equal 
heights, lower than that of i, the particle is allowed to 
move to the top of cither neighbour with equal probabil- 
ity. 

To extend the model to larger diffusion lengths, a parti- 
cle arriving at a chosen site is allowed to diffuse to farther 
neighbours on either side. As an illustration, for allowed 
diffusion length Nd = 2, a particle arriving at a chosen 
site i is allowed to diffuse up to two neighbours, succes- 
sively lower, on either side. The particle first moves to 
the site lowest among i, i + I and i — I. If it has already 
moved to a neighbouring site, i.e., i — 1 or i H- 1, it then 
moves to the next neighbour, i.e., i — 2 or i + 2 respec- 
tively, if that site is further lower. If in the first move 
i — 1 and i + 1 sites are equal to each other and lower 
than site i, the movement to the left or right is initiated 
by random choice. If however both sites i — 1 and i + 1 are 



equal to or higher than site i, no movement takes place. 
Figure [T] illustrates the deposition rules for the model. 



FIG. 1: Schematic diagram for deposition rules in 

RDSD. The numbers (a/b) at the top of columns show 

maximum possible left diffusion (a) and right diffusion 

(b) at that site. 



III. ANALYTICAL ESTIMATE OF 
DIFFUSIBILITY 

Let us consider a system of length L. The probability 
of selecting a random site among the L available sites 
is p = J. Let us first focus on the morphology of a 



FIG. 2: A configuration where at least m step diffusion 
on both sides is possible. 

small subsystem consisting of the chosen site i of height 
h (where a new particle arrives) and its two nearest neigh- 
bours {i — 1) and (z -|- 1) with heights ai and 5i respec- 
tively. The configurations which allow at least one step 
diffusion at the site i, say, to the left, must have the 
height ai of the {i — l)th column less than the height h 
of the ith column, i.e., ai < h. If n layers, or a total 
of nL particles are deposited, the probability that the 
chosen site has a lower left neighbour is given by 



Hi 
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The summations are subject to the restrictions ai < h 
and ft. -I- ai < nL. Likewise, for diffusion to both left 
and right, at the site i, by at least one step, the height 
of the left nearest neighbour ai and that of the right 
nearest neighbour bi should both be less than h. Then 
the probability for at least one step diffusion to both sides 



is (See Figure HJ, 

nL nL—h nL—h—ai 
7 = 1 ai=0 hi=C 
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Here again, the summations are subject to the restric- 
tions ai < h, hi < h and /i + ai + 6i < nL. From the 
set of aU configurations that allow at least one step dif- 
fusion to the left at the chosen site, if one removes all 
those that allow diffusion to both sides by at least one 
step, one is simply left with the configurations that allow 
diffusion only to the left and not to the right. Hence, the 
probability of diffusing at least one step to the left with 
no diffusion to the right is given by 
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Hi - H^^'-^^ 



From symmetry, we may write T\ , the probability of 
diffusion to the right, with no diffusion to the left as, 

r(^) = Hi-H(^'^). 

The probability that a particle arriving at the chosen 
site may diffuse to at least one neighbouring site on either 
side is given by. 



(L) 



Ai = T\^' + V\ 
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Extending the above arguments to m diffusion steps, the 
probability of diffusion by at least m steps to cither side 
is given by, 

A™ = 2H™ - H(,f ^^' m = l,2,... (5) 

where, following Eq. ^ we may write 

\ ^ nL /^ nL — h/^ 

m — / ^ '^h <-'ai 

h,{ak} 

nL—{h+ai^ hOfc^i)^ nL—{h+ai^ |-a„_i)^ 

h+aiA hOmrj^ -^ fjji^ l\p]nL-{h+ai^ ha„) 

where, the summations are over h and every Ofc, for k = 
1, 2, . . . , 771. The restrictions on the summations are < 
o-rn < im-i < • ■ • < fli < /i and h + ai + ■ ■ ■ + a„i < nL. 
Similarly, we may write for H,„ ' , 

h,{ak.bk} 

nL— (/i+aiH hat-i)^ nL— (/i+aiH ha^^i)^ 

'-'a.k ■ ■ ■ '^am 

nh — hf^ nL—{h-\-bi) /^ nL—{h-\-bi-\ \-bni-i) /^ 

/i+aiH ha„+6iH h6,„ 



where, the restrictions on the summations are < Om < 
Um-i < ■ ■ ■ < ai < h, < bm < bm-1 < ■ ■ ■ < bi < h 
and h + ai + ■ ■ ■ + a„i + bi + ■ ■ ■ + bm ^ nL. 

The probability that a particle arriving at a selected 
site will diffuse exactly ni steps is given by 

"Pm = A™ - A™+i m = 1, 2, . . . (6) 

and the probability of zero or no diffusion is 

n = l-Ai. 

We have computed these probabilities for several sys- 
tem sizes. In Table U and Table [H] the values are given 
for system sizes L = 48 and i = 512 respectively. We 
observe that the system sizes do not have any marked 
influence on the probabilities. 

A plot of Vrn versus the number of layers in logarith- 
mic scale is shown in Fig [3] The probability that there 
is no diffusion at all, steadily decreases as the number 
of deposited layers increases. For finite diffusion lengths, 
all other probabilities Vm, ni ~ 1,2,... increase as the 
number of layers increases. After 10 layers of deposi- 
tion, the zero and one step diffusion probabilities are the 
highest, followed by two and three step diffusion prob- 
abilities. After 100 layers, when the RDSD model has 
attained saturation, the nearest neighbour diffusion prob- 
ability dominates, followed by zero diffusion probability 
as can be seen in Fig |4l The probabilities for diffusion 
lengths 2 and 3 are significantly smaller in comparison, 
together accounting for less than 1/3 of the total possi- 
bility of diffusion. The probabilities for higher diffusion 
lengths are negligible. 
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FIG. 3: (Color Online) Probabilities of diffusion to zero, 
one and two steps plotted versus number of layers in log 

scale. 

In any given layer, beyond the first, the diffusion prob- 
ability maximizes at diffusion length one. Thus, one step 
diffusion has the most dominant effect in determining 
surface smoothness. 
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Vo 


0.5253 


0.4533 


0.4035 


0.3810 


0.3672 


0.3478 


Pi 


0.3882 


0.4048 


0.4050 


0.4016 


0.3970 


0.3902 


V2 


0.08207 


0.1252 


0.1579 


0.1728 


0.1824 


0.1918 


Vs 


0.004444 


0.01578 


0.02998 


0.03855 


0.04499 


0.05771 


Vi 


0.00004930 


0.0009140 


0.003366 


0.005471 


0.007388 


- 


Vs 


7.166 X 10"* 


0.00002252 


0.0002457 


0.0005842 


0.0008818 


- 



TABLE I: Probability Vr, 



of a site with maxinium diffusibility m on cither side after n 
particle deposited for system size i = 48. 



1,2,5,10,20, 100 layers of 



No. of Layers 
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5 
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20 


100 
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0.5279 


0.4542 


0.4172 


0.3816 


0.3673 


0.3488 


Pi 


0.3864 


0.4046 


0.4005 


0.4011 


0.3970 


0.3900 


V2 


0.08105 


0.1247 


0.1514 


0.18664 


0.1824 


0.1914 
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0.004516 


0.01566 


0.02762 


0.02452 


0.04497 


0.05747 


Vi 


0.00006098 


0.0009081 


0.002969 


0.005475 


0.007382 


0.01064 


Vs 


1.623 X 10"'^ 


0.00002336 


0.0002037 


0.0005439 


0.0008808 


- 



TABLE IL Probability P„ 



of a site with maximum diffusibility m on either side after n 
particles deposited for system size L = 512. 



1,2, 5, 10, 20, 100 layers of 
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FIG. 4: (Color Online) Diffusion probabilities versus 

diffusion length for various layers. For more than two 

layers, the probabilities maximize for nearest neighbour 

diffusion. 



IV. NUMERICAL SIMULATION STUDY 

Our simulations for nearest neighbour or one step dif- 
fusion include all system sizes studied by Family, L — 
24, . . . ,384 and beyond, up to L = 3072. Limiting our 
study to the system sizes studied by Family in (1 + 1) di- 
mension reproduces his results with great accuracy. The 
study is then extended to larger diffusion lengths. For 
diffusion lengths Nd ~ 2,3,4,6,8--- we observe that a 
small change in the surface properties is discernible up to 
Nd = 3. The surface roughness decreases slightly, system 
saturates at an earlier time and has a lower saturation 
width. Beyond Nd = 3, the interface width versus time 
plots merge (See Figure [S]). 



FIG. 5: (Color Online) \nWit) vs Int with i = 48 for 
different diffusion lengths. 



It is important to note that in our simulations of 
RDSD, two independent random number generators were 
used, one for selecting the site on the growing surface and, 
another for depositing the particle with equal probability 
when left and right diffusions are equally likely. These 
two random number generators arc completely indepen- 
dent and uncorrclated to each other and are generated 
by linear congruential method^ii. 

To understand the results of simulation better, we 
study the maximum possible diffusibility of sites at vari- 
ous stages of growth. The maximum possible diffusibility 
for a given site is the number of adjacent sites which are 
successively lower in height (See Figure [T]). For RD, we 
find that beyond sub monolayer region, about 70% of the 
sites allow diffusion to at least one step or more, implying 
that an arriving particle if allowed, will diffuse to one or 



more, successively lower nearest neighbours (See Figure 
ini . The plots are qualitatively similar to the plots in Fig- 
ure [31 The percentage of sites with a certain diffusibility, 
say m, is a measure of the probability of diifusion by m 
steps, and the number of layers deposited is a measure 
of time. The diffusibility distribution is next studied for 




FIG. 6: (Color Online) Number of sites (in %) with 

different diffusibility versus time for Random 

Deposition Model. 




for first few layers most 

sites have small deviation aboui 

a mean height 



with increase in layers 
more deviations appear 
but the distribution 
saturates 
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FIG. 7: (Color Online) Height distribution in simulation 

of nearest neighbour diffusion model for system size 

L = 48. The distributions are sharp and narrow for 

t < 403 (growth region). The flatter, wider curves for 

t > 403, overlap, showing saturation of the distribution. 



Family's nearest neighbour model where a newly arriv- 
ing particle diffuses to a lowest nearest neighbour, even 
when farther lower sites are available. The height profile 
is studied after saturation, to ascertain the possibility of 
further diffusion at each site. We observe that in the 
saturation region, for about 9% of the sites, further dif- 
fusion, (beyond nearest neighbour) is possible. With the 
diffusion length extended to 2, it is observed that, for 
only 0.8% of the sites, diffusion beyond two neighbours 
is possible. With further increase in Nd, the percentage 
of sites favouring diffusion greater than the maximum 
allowed limit set by Nd, diminishes further. These ob- 
servations agree with the diminishing probabilities Vm 
appearing in Eq. (jS]), for increasing m. The analytically 
calculated probabilities and the diffusibility distributions 
in the simulations explain the small change observed in 
the growth exponent and saturation width when the dif- 
fusion length is increased from 1, 2 to 3 steps. Further 
increase in diffusion length beyond 3 has no discernible ef- 
fect on the surface properties obtained from simulations. 
A study of the height profile, shows that in RD, the dis- 
tribution of the heights of the sites about a mean height, 
becomes flatter and wider as more layers are deposited, 
signifying an increase in interface roughness as time pro- 
gresses. When diffusion is allowed, this flattening and 
widening process stops beyond a certain time, implying 
saturation of interface width. In Figure [71 this saturation 
of the height distribution is shown for diffusion length 
Nd = 1 and system size L = 48. The height distributions 
after saturation in RDSD models depend on the diffusion 
length and arc plotted for Nd = 1, 2, 3, 4, 8 and 16 in Fig- 
ure |H1 It is noted that the saturated height distributions 



do not change for diffusion length Nd beyond 3. 
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"Saturation height 
distribution" changes 
with diffusion, but with 
further increase in 
diffusion steps no 
further change is seen 
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FIG. 8: (Color Online) Height distribution after 

saturation in simulations of diffusion lengths = 2, 3, 4, 8 

and 16 for system size L = 48. 

In RDSD, the logarithmic graph of interface width 
versus time, is initially, nearly coUinear with the corre- 
sponding graph in RD; later it deviates from the lin- 
ear form and eventually saturates. We estimate the 
quantitative nature of this deviation. Let us consider 



a random deposition process with system size L, num- 
ber of particles deposited N and N < L defines sub- 
monolayer deposition. Tlic probability of a site being 



occupied is -^ and that being empty is 1 — -j-. In the 
sub-monolaycr random deposition, a chosen site may 
allow diffusion provided either one or both its nearest 
neighbours are empty. Thus, the probability of find- 
ing a site where diffusion may be possible is given by 

(f)(f)(l-f) + (f)(l-f)(f) + (f)(l-f)(l-f) 
which is equal to f (l - f )(l + f ). 



ior t < -j= and the deviation will become significantly 

large for t > -y=. 

In Figure ini the results of an actual simulation for L = 
48, Nd = 1,2 and RD are shown where departure from 
random behaviour is evident for In t > In ^ . 

v3 



DISCUSSION AND CONCLUSION 
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FIG. 9: (Color Online) InlF(t) vs Int showing deviation 
from the straight line nature. 

The effects of diffusion are most noticeable when this 
probability function is maximum. Denoting ^ = t, it is 
easily verified that the maximum for the above probabil- 
ity function f {t) = t (1 — P) occurs when t ~ -y^- Thus 

for t ~ -i=, fmax — 0.4, which implies no more than 40% 
of the sites in the sub-monolayer in RD, allow diffusion. 
In an actual process, which involves diffusion, we conjec- 
ture that the above probability function will reach the 
maximum value at an even lower value of t. Deviation 
from random behaviour, in an actual process, is visible 



We have shown by means of analytical calculation as 
well as numerical simulation that in RD, appearance of 
long staircase-like formations with successive sites having 
ascending or descending heights, even after several lay- 
ers of deposition, is improbable. After a few successive 
descending sites (1 or 2), one reaches the local minimum 
of a staircase. Thus, even in the case of RD, maximum 
diffusibility of sites is limited. However, the differences in 
the heights of unequal nearest neighbours increase with- 
out limit. This is in agreement with the fact that rough- 
ness of the interface increases as i^/^, even though the 
maximum diffusibility at any site is liiTiited. It is in- 
teresting to note that the number of sites with different 
diffusibility approaches stable values as time progresses 
(See Figures 13] and [ni . 

The study of the effect of diffusion can be extended 
to systems that allow diffusion through flat regions or 
even through barriers. They may represent continuation 
of rolling of classical particles through flat regions due 
to inertia or quantum barrier penetration, e.g., chemical 
deposition process even for reactants separated by po- 
tential barrier. Further investigation pertaining to the 
above processes are in progress. 
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